Impact of string and monopole-type junctions on domain wall dynamics: implications 

for dark energy 
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We investigate the potential role of string and monopole-type junctions in the frustration of 
domain wall networks using a velocity-dependent one-scale model for the characteristic velocity, v, 
and the characteristic length, L, of the network. We show that, except for very special network 
configurations, v 2 < (HL) 2 < (p a + p M )/p m where H is the Hubble parameter and p CT , p M and 
pm are the average density of domain walls, strings and monopole-type junctions. We further 
show that if domain walls are to provide a significant contribution to the dark energy without 
generating exceedingly large CMB temperature fluctuations then, at the present time, the network 

2/3 K 2/3 

must have a characteristic length Lo < 10f2 CTO kpc and a characteristic velocity vo < 10 ^ CT o 
where fl CT o = pao/pco and p c is the critical density. In order to satisfy these constraints with 
^o-o ~ 1, PmO would have to be at least 10 orders of magnitude larger than p a o, which would be in 
complete disagreement with observations. This result provides very strong additional support for 
the conjecture that no natural frustration mechanism, which could lead to a significant contribution 
of domain walls to the dark energy budget, exists. 



I. INTRODUCTION 

There is now very strong observational evidence that 
our Universe is currently undergoing a period of acceler- 
ated expansion [3, 0] . In the framework of General Rel- 
ativity, the observed acceleration of the Universe ought 
to be explained by the existence of an exotic dark energy 
component which violates the strong energy condition. 
In [J] it was suggested, for the first time, that a frozen 
domain wall network could be responsible for such ac- 
celeration and several other authors have subsequently 
advocated this possibility 0, [ESQ- 

The evolution of cosmological domain wall networks 
has been studied in detail using both high-resolution 
numerical simulations and a semi-analitical velocity- 
dependent one-scale (VOS) model @, ©, EE El El El • 
Strong analytical and numerical evidence was provided 
that a frustrated domain wall network, accounting for a 
significant fraction of the energy density of the Universe 
at the present time, will never emerge from realistic phase 
transitions. These results appear to rule out a significant 
contribution of domain walls to the dark energy budget. 

In all these studies the energy of the domain walls was 
assumed to be very small. Still, domain wall junctions 
can be responsible for freezing a domain wall network if 
they are heavy enough. In that case their contribution 
to the energy budget could not be neglected but it has 
been argued in [l(J E3] that it would spoil the dark en- 
ergy properties associated with a frustrated domain wall 
network. More recently, a mechanism for dynamical frus- 
tration of domain wall networks involving kinky vortons 
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was suggested in [14j. 

In this paper we investigate the impact of string and 
monopole-type junctions on the dynamics of domain wall 
networks. Accounting for the detailed contribution of the 
junctions on domain wall network simulations is not a 
trivial task. In fact, the standard Press- Ryden-Spergel 
(PRS) algorithm [15|] (implemented in all field theory 
numerical simulations in order to ensure fixed comoving 
resolution) increases artificially the impact of the junc- 
tions on the overall network dynamics during the course 
of the simulations. This effect is not very important for 
the light junctions which are usually considered in such 
simulations but could be relevant in the case of heavy 
junctions. 

In order to overcome this potential problem we develop 
a semi-analytical VOS model which incorporates the con- 
tribution of string and monopole-type junctions to the 
overall dynamics of domain wall networks. We use this 
model to constrain the characteristic length and velocity 
of domain wall networks and discuss the corresponding 
cosmological implications, in particular for dark energy. 



II. SEMI-ANALITICAL VOS MODEL WITH 
JUNCTIONS 

Consider the configuration presented in figure [1] where 
4 straight string-type junctions of energy per unit length 
H (represented by 4 large dots) are connected by 4 pla- 
nar domain walls of energy per unit length a (repre- 
sented by 4 straight lines) defining a square domain of 
characteristic size I (Jo represents the characteristic size 
prior to collapse). Here we assume that nothing varies 
along the direction perpendicular to the square domain, 
so that the infinite string-type junctions have no curva- 
ture. In this case the domain wall dynamics is effectively 
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FIG. 1: The left panel represents the domain wall configu- 
ration prior to collapse. Solid lines represent domain walls 
of superficial tension a and the large dots represent string 
junctions of tension p. The right panel represents the config- 
uration after some time has elapsed. 



2-dimensional. The energy per unit length of this config- 
uration is given by 



E = iajl + 4/i 7At + 2V2o-(/ - , 



(1) 



where 7 = 7(r = (1 - v 2 )" 1 ' 2 , 7/i = (1 - vD' 1 ' 2 , v 
v a = — {dl/dt)/2 represents the domain wall velocity and 
= \/2v is the velocity of the string-type junctions. 
For simplicity, we shall make the approximation that 
the shape is maintained during the collapse of the square 
domain so that the curvature of the domain walls remains 
concentrated at the vertices. Although we expect domain 
walls to acquire curvature in the course of the collapse, 
this approximation will capture the relevant physics and 
does not affect our main results. It follows from energy 
conservation that, in Minkowski spacetime, the equation 
of motion is given by 



dv 
~db 



(1 
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and 
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(2) 



(3) 



(4) 



These results do not directly apply to the case of a do- 
main wall network since they were obtained for a very 
particular configuration. Still, we shall use Eq. |2|) as a 
guide into the fundamental aspects of network evolution. 

The characteristic lengths, and L a , of a two- 
dimensional domain wall network with junctions are de- 
fined as 
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(5) 



where p M and p a represent the average density of string- 
type junctions and the average domain wall density, re- 
spectively. The relation between L M and L a depends on 
the geometrical properties of the domain wall network. 
For a regular hexagonal lattice with Y-type junctions 
L^/L a = 3 1 / 4 , for a regular square lattice with X-type 
junctions L^/L^ = 2, while L^/ L a = 2 1 / 2 3 3 / 4 for a reg- 
ular triangular lattice with *-type junctions. In general, 
L^/La ~ 1 and, consequently, in the following we shall 
assume that L = = L a . The characteristic velocities 
are also expected to be similar and consequently we shall 
make the approximation that v = v^ = v a . 

Then substituting I = L/2 in Eq. ([2]) one obtains, in 
the non-relativistic limit, 



dv 
~dt 



1 



k 



= 0. 



(6) 



with k = 2s (0) = 1 — \/2/2. Note that the domain wall 
configuration presented in Fig. 1 is very specific. In 
general the curvature parameter, k, depends on the geo- 
metrical properties of the domain wall network [l2j • 
Eq © can also be re-written as 



dv 



L 



(1-/m) = 0, 



(7) 



where / M = p^/p and p = p a + p^. 

The generalization to the 3-dimensional case is triv- 
ial. In this case we need to consider the impact of 
monopole-type junctions, of mass m and average density 
p m = m/L^, on the network dynamics. The evolution 
equation for the characteristic velocity becomes, up to 
coefficients of order unity, 



dv k , „ . 
~dt ~L = 



(8) 



The 



with L = L m = = L a and v = v m — — 
energy density fractions of the various components are 
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(9) 
(10) 

(11) 



where p = p a + p M 

We have implicitly assumed that the energy of the 
domain walls as well as that of string and monopole 
type junctions is localized, which may not always be a 
good approximation. For example, in the case of global 
monopoles described by a Lagrangian with a standard ki- 
netic term there is a nearly constant long-range force be- 
tween monopole and anti-monopole pairs which could be 
dynamically relevant. However, such long-range forces, 
if they are important, would trigger a very effective an- 
nihilation between monopole and anti-monopole pairs. 
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This would constitute an additional obstacle to the frus- 
tration of a domain wall network and consequently we 
shall not consider them in the present paper. Long range 
forces between monopoles can be eliminated by consider- 
ing non-standard kinetic terms which localize the energy 
of the monopoles inside their core. 

In an expanding universe we also need to account for 
the deceleration caused by the Hubble expansion. In the 
case of planar domain walls, the momentum per comov- 
ing area is proportional to a" 1 , so that 

la v„(xa- 3 &^+3v a H(l-vl)=0. (12) 

Similarly, for cosmic strings the momentum per comoving 
length is proportional to a -1 . Consequently 

7 ^ M cxa- 2 o^-+2tyff(l-i£)=0. (13) 

For massive junctions one has 

j m v m (xa- 1 ^^+ v m H(l - v 2 m ) = . (14) 

In the absence of energy loss mechanisms the evolution 
of the average energy density is given by 

Pa^^a" 1 ^+fr(l + 3t£K=0, (15) 

/V« 7 ^ 2 ^+2(1 + 1^)^ = 0, (16) 

p ro oc 7ro a- 3 «■ ^ + 3(l+^\ Hp m = 0,(17) 



the non-relativistic limit this term is negligible. 

In any case, it is easy to verify that this term leads to a 
larger L and, consequently, it does not help frustration. 

We use the energy density fractions as weight factors 
in the estimate of the Hubble damping term which needs 
to be added to Eq. Q . The equation describing the evo- 
lution of the characteristic velocity of the network then 
becomes 



0, (20) 



Solving Eq. (|20p in the matter dominated era, assum- 
ing that = and p m p a and neglecting the decaying 
mode, one obtains 



^ 2 dv 2 k p a 
5 dt 7 L p m 



(21) 



A similar analysis for p a = and p m 3> p^ would give 

Hv= 2dv = 2kp IL 

3 dt 5Lp m y ' 

In both limits one has 



v 2 <HLv~k Pa + P » , (23) 



where we have taken into account that v < LH since the 
characteristic velocity does not change abruptly. 

III. IMPLICATIONS FOR DARK ENERGY 



for domain walls, cosmic strings and point masses. In the 
course of network evolution energy will be exchanged be- 
tween the different components, leading to similar char- 
acteristic velocities (v m ~ ~ v a ). To account for that 
one needs to include extra terms Q a , and Q m on the 
right hand side of Eqs. (|15H17|) , respectively. Here, Qx 
represents the energy transfered per unit of time and vol- 
ume from the component X to the other two components 
(note that Q m + + Q a = 0). 

The evolution of the total energy density is given by 

§+[(! + U + 2/-) + « 2 (3 - U - 2/*)] Hp = , (18) 



where we have used Eqs. (|15H17p , taking into account 
that p = p a + p^ + p m and the fact that the Qx terms 
cancel out. 
Hence, 



dL 
~dt 



A + U + 2U 
3 - /„ - 2f a 



HL, 



(19) 



so that for v = v a = = v m <C 1 one has dL/dt = HL 
(L oc a) with L = L a = L „ = L m . A term proportional 
to vp/L could be added to the right hand side of Eq. (TT8")) . 
in order to account for energy losses by the network. In 



The amplitude of the fractional energy density fluc- 
tuations, 5, associated with domain walls, on a physical 
scale, Ly, much larger than the characteristic scale, L, 
of a domain wall network, is given approximately by 



Spa 

Pc 



<1n 



N 



(24) 



where N ~ (Ly/L) 3 is the number of domain walls on a 
volume V = Ly, fl a = p a /p Cl p c is the critical density 
and 8 p a is the root mean square fluctuation of the do- 
main wall energy density on the physical scale Ly . The 
amplitude of CMB temperature fluctuations generated 
by domain walls, around the present time, is constrained 
to be smaller than 10 -5 down to scales of the order of 
Ly = -ff,7 1 /100 (in order not to spoil the agreement be- 
tween the theoretical model and observations 0,13). This 
implies 



10 3 n a0 {H Q L f 2 < 10" 



(25) 



Consequently, HqLq < 10 5 f^Q ^ 3 which results in the 
2/3 

constraint Lq < 10 fl aQ kpc. Unless there is an abrupt 
velocity change, this also translates into a stringent limit 
on the characteristic velocity of the domain walls at the 
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present time, vq < H L < 10~° £l a0 ' . Hence, using 
Eq. (f2"5)) one obtains 

k PaO±P^ < (i/oLo)2 < 1Q _ 10 Q -4/3 (26) 

The value of the curvature parameter has been estimated 
using high-resolution numerical simulations of domain 
wall evolution. For standard domain wall networks with- 
out junctions k ~ 1 while for domain wall networks with 
junctions a smaller value has been observed, but still of 
order unity [l2| ■ A value of k <C 1 is only expected in the 
case of very special configurations such as 2-dimensional 
hexagonal lattices with Y-type junctions, square lattices 
with X-type junctions and triangular lattices with *-type 
junctions. These are very unnatural configurations cor- 
responding to very specific initial conditions which would 
violate causality if they were to extend over scales larger 
than the particle horizon. Consequently we shall assume 
that k ~ 1. 

Note that domain walls would need to have an aver- 
age energy density of the order of the critical density to 
provide a significant contribution to the dark energy. We 
conclude, from Eq. (|26p . that the average energy density 
of the monopole-type junctions would have to be 10 or- 
ders of magnitude greater than that. Such a high value 
of the energy density is in complete disagreement with all 
cosmological evidence. Of course, if £l a o was very small 
then frustration could effectively occur but, in that case, 
domain walls would not play any relevant role as a dark 
energy component. 

IV. CONCLUSIONS 

In this paper we have imposed strong constraints on 
the characteristic length and velocity of domain wall net- 



works with string and monopole-type junctions using a 
semi-analitical VOS model. We have shown that a suc- 
cessful domain wall scenario for dark energy would re- 
quire that Lo < lOkpc and vq < 10~ 5 . We have demon- 
strated that such small values of Lq and v o could only be 
achieved if the contribution of the monopole-type junc- 
tions to the total density of the universe was several or- 
ders of magnitude larger than that of domain walls and 
strings (assuming fl a Q ~ 1), in complete disagreement 
with observations. 

These results highlight the main difficulty with alterna- 
tive mechanisms for the frustration of domain wall net- 
works. The inclusion of additional degrees of freedom 
such as heavy junctions and friction may slightly reduce 
the characteristic length and velocity of the domain walls 
but is insufficient to lead to frustration, due to the lim- 
ited amount of matter with which domain walls can in- 
teract while conserving energy and momentum at the 
present time (the mechanism for the frustration of do- 
main wall networks proposed in [14J is expected to face 
similar problems). Our results provide further evidence 
for the absence of a successful mechanism which can lead 
to the frustration of cosmologically relevant domain wall 
networks. 
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